STABLE COMPLETE INTERSECTIONS 



LORENZO ROBBIANO AND MARIA LAURA TORRENTE 

Abstract. A complete intersection of n polynomials in n indeterminates has 
only a finite number of zeros. In this paper we address the following ques- 
tion: how do the zeros change when the coefficients of the polynomials are 
perturbed? In the first part we show how to construct semi-algebraic sets 
in the parameter space over which all the complete intersection ideals share 
the same number of isolated real zeros. In the second part we show how to 
modify the complete intersection and get a new one which generates the same 
ideal but whose real zeros are more stable with respect to perturbations of the 
coefficients. 



1. Introduction 

What is the defining (or vanishing) ideal of a finite set X of points in the affine 
space? The standard answer is that it is the set of all the polynomials which vanish 
at X. And there are very efficient methods to compute it, based on Buchberger- 
MoUer's algorithm (see for instance [1], [2] and [6]). 

However, the logical and computational environment changes completely when 
the coordinates of the points are perturbed by errors, a situation which is normal 
when dealing with real world problems. In that case one has to use approximation 
and to consider the question of stability. Introductory material about this topic 
can be found in the book [3 , in particular in the paper |14j and its bibliography. 

The methods used so far share the strategy of modifying the Buchberger-MoUer 
Algorithm and compute a Grobner basis or a border basis of an ideal of polyno- 
mials which almost vanish at X (see for instance flOj and [11]). A key remark is 
that, whatever algorithm is used, at a certain moment one has computed n poly- 
nomials /i, . . . ,/ti which generate a zero-dimensional ideal. Since the dimension 
has dropped from n to zero, the n polynomials form a complete intersection which 
almost vanishes at X. Further steps in the algorithm will be used to eliminate 
spurious points and to produce a Grobner or border basis. 

Now, a complete intersection of n polynomials in n indeterminates has only a 
finite number of zeros, and the main question is: how do the zeros change when the 
coefficients of the polynomials are perturbed? Can we devise a strategy to make the 
situation reasonably stable? In other words, can we change the generating polyno- 
mials so that the stability of their common zeros increases? It is well-known that for 
a linear system with n equations and n unknowns, the most stable situation occurs 
when the coefficient matrix is orthonormal. Is there an analogue to orthonormality 
when we deal with polynomial systems? 
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In numerical analysis the condition number of a problem measures the sensitivity 
of the solution to small changes in the input data, and so it reveals how numerically 
well-conditioned the problem is. There exist a huge body of results about condition 
numbers for various numerical problems, for instance the solution of a linear system, 
the problem of matrix inversion, the least squares problem, and the computation 
of eigenvalues and eigenvectors. 

On the other hand, not very much is known about condition numbers of polyno- 
mial systems. As a notable exception we mention the paper [17j of Shub and Smale 
who treated the case of zero-dimensional homogeneous polynomial systems; later 
on their result was extended by Degot (see [8 ) to the case of positive-dimensional 
homogeneous polynomial systems. 

Tackling the above mentioned problem entails a preliminary analysis of the fol- 
lowing question of algebraic nature. If we are given a zero-dimensional complete 
intersection of polynomials with simple zeros, how far can we perturb the coeffi- 
cients so that the zeros remain smooth and their number does not change? It is 
quite clear that smoothness and constancy of the number of zeros are essential if 
we want to consider the perturbation a good one. 

Starting with the classical idea that a perturbed system is a member of a family of 
systems, we describe a good subset of the parameter space over which the members 
of the family share the property that their zero sets have the same number of 
smooth real points. This is the content of Sectio n [2| w here we describe a free 
(see Proposition 2.6), and a smooth (see Theorem 2.12) locus in the parameter 
space. Then we provide a suitable algorithm to compute what we call an /-optimal 
subscheme of the parameter space (see Corollary |2.16[ ): it is a subscheme over 
which the complete intersection schemes are smooth and have the same number of 
complex points. The last important result of Section[2]is Theorem 2.20 which proves 
the existence of an open non-empty semi-algebraic subscheme of the /-optimal 
subscheme over which the number of real zeros is constant. 

Having described a good subscheme of the parameter space over which we are 
allowed to move, and hence over which we can perturb our data, we pass in Section[3] 
to the next problem and concentrate our investigation on a single point of the zero 
set. After some preparatory results, we introduce a local condition number (see 
Definition 3.14) and with its help we prove Theorem 3.15 which has the merit of 
fully generalizing a classical result in numerical linear algebra (see Remark 3.161. 

The subsequent short Section [4] illustrates how to manipulate the equations in 
order to lower, and sometimes to minimize, the local condition number (see Propo- 
sition 4.1 1. Then we concentrate on the case of the matrix 2- norm and show how 



to achieve the minimum when the polynomials involved have equal degree (see 
Proposition 4.3). The final Section [5] describes examples which indicate that our 
approach is good, in particular we see that when the local condition number is 
lowered, indeed the corresponding solution is more stable. 

This paper reports on the first part of a wider investigation. Another paper is 
already planned to describe how to deal with global condition numbers and how to 
generalize our method to the case where the polynomials involved have arbitrary 
degrees. 

All the supporting computations were performed with CoChA (see ) . We thank 



Marie-Frangoise Roy and Saugata Basu for some help in the proof of Theorem 2.20 
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2. Families of Complete Intersections 

Given a zero-dimensional smooth complete intersection X, we want to embed it 
into a family of zero-dimensional complete intersections and study when and how X 
can move inside the family. In particular, we study the locus of the parameter-space 
over which the fibers are smooth with the same number of points as X, and we give 
special emphasis to the case of real points. 

We start the section by recalling some definitions. The notation is borrowed 
from [15] and [16], in particular we let xi, . . . ,a;„ be indeterminates and let T" be 
the monoid of the power products in the symbols Xi, . . . ,a;„. Most of the times, 
for simplicity we use the notation If if is a field, the multivariate 

polynomial ring K[y:] = K[xi, . . . , Xn] is denoted by P, and if /i(x), . . . , /^(x) are 
polynomials in P, the set {/i(x), . . . , /fc(x)} is denoted by f (x) (or simply by f). 
Finally, we denote the polynomial system associated to f (x) by f{x) = (or simply 
by f = 0), and we say that the system is zero-dimensional if the ideal generated 
by f(x) is zero-dimensional (see [TS], Section 3.7). 

Easy examples show that, unlike the homogeneous case, in the inhomogeneous 
case regular sequences are not independent of the order of their entries. For in- 
stance, if /i = y{x + 1), /2 = z{x + 1), /s = X, then (/i,/2,/3) is not a regular 
sequence, while (/3,/i,/2) is such. However, we prefer to avoid a distinction be- 
tween these cases, and we call them complete intersections. In other words, we use 
the following definition. 

Definition 2.1. Let i be a positive integer, let f(x) be a set of t polynomials in 
P = K[xi, . . . , Xn] and let / be the ideal generated by f (x). 

(a) The set f(x) (and the ideal /) is called a complete intersection if the 
equality dim(P//) = n ~ t holds. 

(b) The set f(x) (and the ideal /) is called a zero-dimensional complete 
intersection if it is a complete intersection and t — n. 

Let n be a positive integer, let P denote the polynomial ring K[xi, . . . ,Xn], let 
f(x) = {/i(x), . . . , /„(x)} be a zero-dimensional complete intersection, and let / 
be the ideal of P generated by f(x). We let m be a positive integer and let 
a = (oi , . . . , am) be an m-tuple of indeterminates which will play the role of param- 
eters. If Fi(a, x), . . . , Fri(a, x) are polynomials in if [a, x] we let i^(a, x) = be the 
corresponding family of systems of equations parametrized by a, and the ideal gener- 
ated by i^(a, x) in if [a, x] is denoted by i(a, x). If the scheme of the a-parameters 
is S, then there is a if-algebra homomorphism : if [a] — > A'[a, x]//(a, x) or, 
equivalently, a morphism of schemes $ : Spec(if [a, x]/i(a, x)) — > S. 

Although it is not strictly necessary for the theory, for our applications it suffices 
to consider independent parameters. Here is the formal definition. 

Definition 2.2. If 5 = and i(a, x) n if [a] = (0), then the parameters a are 
said to be independent with respect to F(a, x), or simply independent if the 
context is clear. 

The first important step is to embed the system f (x) = into a family, but we 
must be careful and exclude families of the following type. 

Example 2.3. Consider the family F(a,x) = {xi{ax2 + l),a;2(aa;2 + !)}• It is a 
zero dimensional complete intersection only for a = while the generic member is 
positive-dimensional. 
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Definition 2.4. Let f (x) be a set of polynomials in P = K[xi, . . . , Xn] so that f (x) 
is a zero-dimensional complete intersection and let F(a, x) be a family parametrized 
by m independent parameters a. We say that i^(a, x) (and similarly if [a, x]//(a, x) 
and Spec(iir[a, x]//(a, x))) is a generically zero-dimensional family contain- 
ing f(x), if f(x) is a member of the family and the generic member of the family 
is a zero-dimensional complete intersection. 

A theorem called generic flatness (see j9j, Theorem 14.4) prescribes the ex- 
istence of a non-empty Zariski-open subscheme U of S over which the morphism 
^^^{U) — > U is flat. In particular, it is possible to explicitly compute a subscheme 
over which the morphism is free. To do this, Grobner bases reveal themselves as a 
fundamental tool. 

Definition 2.5. Let F(a, x) be a generically zero-dimensional family which con- 
tains a zero-dimensional complete intersection f(x). Let S = be the scheme 
of the independent a-parameters and let $ : Spec(iir[a, x]//(a, x)) — > S be the 
associated morphism of schemes. A dense Zariski-open subscheme hi oi S such 
that ^^^iU) — >U is free (flat, faithfully flat), is said to be an /-free (/-flat, 
/—faithfully flat) subscheme of S or simply an /-free (/-flat, /-faithfully flat) 
scheme. 

Proposition 2.6. With the above assumptions and notation, letl{a,x) be the ideal 
generated by /"(a, x) in /^[a, x], let a be a term ordering on T" , let G'(a, x) be the 
reduced a-Grdbner basis of the ideal /(a, x)/C(a)[x], let d{a) be the least common 
multiple of all the denominators of the coefficients of the polynomials in G(a, x), 
and let T = T" \ LT^(/(a, x)/f (a)[x]). 

(a) The open subscheme 14 of defined by d{a) ^ Q is I -free. 

(b) The multiplicity of each fiber over lA coincides with the cardinality of T . 

Proof. The assumption that /"(a, x) is a generically zero-dimensional family im- 
plies that Spec(/i'(a)[x]//(a, x)/ir(a)[x]) — > Spec(/ir(a)) is finite, in other words 
that /C(a)[x]//(a, x)/i'(a)[x] is a finite-dimensional /C(a)-vector space. A standard 
result in Grobner basis theory (see for instance [15j, Theorem 1.5.7) shows that the 
residue classes of the elements in T form a /f'(a)-basis of this vector space. We 
denote by U the open subscheme of A^ defined by d{a) ^ 0. For every point in 
the given reduced Grobner basis evaluates to the reduced Grobner basis of the cor- 
responding ideal. Therefore the leading term ideal is the same for all these fibers, 
and so is its complement T. If we denote by /C[a](j(a) the localization of K[a\ at the 
element d(a) and by /(a, x)"^ the extension of the ideal /(a, x) to the ring K[si\ii{a), 
then /ir[a](i(a) [x]//(a, x)*^ turns out to be a free /ir[a]^(a)-module. So claim (a) is 
proved. Claim (b) follows immediately from (a). □ 

Remark 2.7. We collect here a few remarks about this proposition. First of all, 
we observe that the term ordering a can be chosen arbitrarily. Secondly, for every 
a £ U lei La be the leading term ideal of the corresponding ideal /„. If cr is a 
degree-compatible term ordering, then La is is also the leading term ideal of the 
homogenization /^°™ of /q, (see |in|. Proposition 5.6.3 and its proof). 

Example 2.8. We consider the ideal / = {fi,g) of K[x,y] where fi = x^ — y, 
g = x{x - l){x + l){x - 2){x + 2){x - 3){x + 3){x + I3){x'^ + X + 1). We check that 
I = (A, /2) where = xy^ + 504x^y - ISSxy"^ + Uy^ - 504x2 + 650xy ~ U7y^ - 
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468x + 133?/. It is a zero-dimensional complete intersection and we embed it into 
the family /(a, x) ~ (ax^ — y,g). If we pick a = Lex with y > x and perform the 
computation as suggested by the proposition, we get the freeness of the family for 
all a. Instead, we get the freeness of the family /(a, x) = (oa;"^ — y,f2) for a 7^ 



(see a further discussion in Example 2.141 



Example 2.9. We let P = C[x], the univariate polynomial ring, and embed the 
ideal / generated by the following polynomial x^ — 3x + 2 into the generically zero- 
dimensional family F{a., x) = {aix^ — a2X + a^}. Such family is given by the canon- 
ical if -algebra homomorphism 

:C[a] — >C[a, x]/{ai,a2,a3)/{aix'^ - a2X + a^) 

It is a zero dimensional complete intersection for 

{a e I ai ^ 0} U {a € I ai = 0, 03 ^ 0}. 
It represents two distinct smooth points for 

{a e I ai 7^ 0, a| - Aaias ^ 0}. 
It represents a smooth point for {a e C"^ | ai = 0, a2 7^ 0}. 
It is not a zero-dimensional complete intersection for {at g C'^ | ai = 0, a2 — 0}. 

This kind of examples motivates the following definition. 

Definition 2.10. Let i^(a, x) be a generically zero-dimensional family contain- 
ing a zero-dimensional complete intersection f(x). Let S — be the scheme 
of the independent a-parameters and let $: Spec(iir[a, x]//(a, x)) — > S be the 
associated morphism of schemes. A dense Zariski-open subscheme lA oi S such 
that — > U is smooth, i.e. all the fibers of $~^(Z//) — > U are zero- 
dimensional smooth complete intersections, is said to be an /-smooth subscheme 
of S or simply an /-smooth scheme. 



For instance in Example 2.9 we have the equality S — and the following open 
set U = {a. ^ \ ai ^ 0, — 4aia3 ^ 0} is /-smooth. 

Remark 2.11. We observe that a dense /-smooth scheme may not exist. It suffices 
to consider the ideal / = (a; — 1)^ embedded into the family (a; — a)^. In any event, 
a practical way to find one, if there is one, is via Jacobians, as we are going to show. 

Theorem 2.12. Let F{a.,x) be a generically zero-dimensional family containing a 
zero- dimensional complete intersection f(x). We let S — A^ be the scheme of the 
independent a-parameters, let /(a, x) be the ideal generated by F(a,x.) in K[a,x], 
let /)(a, x) = det(JacF(a, x)) be the determinant of the Jacobian matrix o/F(a, x) 
with respect to the indeterminates x, let J(a, x) be the ideal sum /(a, x) -I- {D{a, x)) 
in K[a^x\, and let H be the ideal in K[a] defined by the equality H = J(a, x)n/i'[a]. 

(a) There exists an I -smooth subscheme of S if and only if H ^ (0). 

(b) If ^ h[a) G H then the open subscheme of S defined by the inequality 
h{a) ^ is I-smooth. 

Proof. To prove one implication of claim (a), and simultaneously claim (b), we 
assume that H ^ (0) and let 7^ h{a) G H . We have an equality of type 
h{a) — a(a, x)/(a, x) -I- &(a, x)/)(a, x) with /(a, x) e /(a, x), and hence an equal- 
^ ~ °ft.^a)'' /(a;x) -I- ^^^^/?(a, x) in J(a, x)/C(a)[x]. For every a G 5 such 
that h{a) ^ the equality implies that the corresponding complete intersection 
has no common zeros with the determinant of its Jacobian matrix, hence it is 
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smooth. Conversely, assume that H = (0). Then the canonical if-algebra homo- 
morphism K[a] — > K[a.,x]/ J{a.,x) is injective and hence it induces a morphism 
Spec (ii'[a, x]/ J(a, x)) — !• of affine schemes which is dominant. It means that 
for a generic point of A^, the scheme Spec {K[a, x]/ J(a, x)) is not empty, and hence 
the corresponding complete intersection Spec (iir[a, x]//(a, x)) is not smooth. □ 

The following example illustrates these results. 

Example 2.13. Let us consider the polynomials fi — xf + x\ - 1, ji 
in C[a;i,a;2] and the ideal I={fi,f2) generated by them. It is a zero-dimensional 
complete intersection and we embed it into /(a, x) = {x^ + aix\ — 1, x\-\- 02X1). 
It is a free family over A^h, and the multiplicity of each fiber is 4. We compute 
D(a, x) = det(JacF(a, x)) and get _D(a, x) = — 2aia2X2 +4x1X2. We let 

J(a, x) = /(a, x) + (_D(a, x)) = (x^ + aiXj — 1, X2 + a2Xi, --2aia2X2 + 4xiX2) 

A computation with CoGcA of Elim([xi, X2], J) yields (^0^02 + 202), and hence 
J(a, x) n K[si\ — {\a\a\ + 2a2). According to the theorem, if U is the comple- 
ment in A^ of the curve defined by \a\a2 + 2a2 = 0, then lA is an /-smooth 
subscheme of A^.. On the other hand, the curve has three components, 02 = 0, and 
aia2 ± 2i = 0. If 02 = then the corresponding ideal is {x\ — 1,X2) which is not 
smooth. If aia2±2i = 0, then the corresponding ideals are (xi=f|^X2 — 1, X2-l-a2Xi) 
which can be written as ((xi ± i)^, x^ -I- 022^1) and hence are not smooth. 

Let us now consider the zero-dimensional complete intersection described by 
the ideal I — (/i, /2) where /i ~ x\ + x\, /2 We embed it into the 

family /(a, x) = (xf — aiXj, x\ + 02X1). As before, it is a free family over A^., 
and the multiplicity of each fiber is 4. We compute D(a, x) = det(jaci;-(a, x)) and 
get i?(a, x) = 2aia2X2 -l-4xiX2. The computation of Elim([x, y], J) yields (0), and 
hence there is no subscheme of A^ which is /-smooth. Indeed, for 02 ^ we have 
/(a, x) = (xi -I- -^x^^ ■^X2 — aix|) which is not smooth. Incidentally, we observe 
that also for a2 — the corresponding zero-dimensional complete intersection is 
not smooth. 

The following example illustrates other subtleties related to the theorem. 
Example 2.14. (Example |2.8| continued) 



We consider the family /(a, x) = [ax^ — y, f2) for a 7^ of Example 2.8 compute 
D{si,yi) = det(JacF(a,x)) and get D{si,yi) = Oax^y^ + 1512ax'* - lOQSax^y + 
126ax2y2 + I950ax3 - SS2ax'^y + y^ + 399ax2 -f lOOSxj/ - ISSy^ - I008x + 650y - 468. 
We let J(a,x) = /(a,x) -I- (/?(a,x)) and get J(a,x) n K[a\ = (/i(a)) where 

9 738170716516748 g 218039463835944563500746 7 166557011563009981474061668 g 
" 7749152384519 " 91409877182005574647 " 31353587873427912103921 " 

276169260891419750846552207 5 986809115998719019081678896 4 63247607413926237871517952 3 
31353587873427912103921 " ^ 31353587873427912103921 " 31353587873427912103921 " 

1316764479863922379654192128 , 317872550804296477704192 974975584016793600000 



31353587873427912103921 13058553883143653521 266501099655992929 

Therefore, if U denotes the complement in K\r of the zeros of /i(a), the theorem 
says that it is a Zariski-open /-smooth subscheme. However, we have already 



seen in Example 2.8 that a = (the origin is va.U) is not in the free locus: we 
observe that the corresponding complete intersection is smooth, but it has only 
two points. The other subtlety is that the Bezout number of the family is 3 x 
4 = 12, but if we substitute y ~ ax^ into /2 we get a univariate polynomial 
of degree 10. The two missing points arc at infinity. No member of the family 
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represents twelve points. The final remark is that if we move the parameter a in 
the locus described by a-h{a) 7^ we always get a smooth complete intersection of 10 
points. If if = C the ten points have complex coordinates, some of them are real, 
but there are no values of a for which all the 10 points are real. The reason is that 
if J'l = — r2 = ~l-^/3» g^j.g ^Yie two complex roots oi + x + l = 0, then two 



of the ten points are {ri,r\), (r2,r2) which are not real points (see Theorem 2.20 



and Example 2.22 1. 



Combining Theorem 2.12 and Proposition 2.6 we get a method to select a Zariski- 
open subscheme of the parameter space over which all the fibers are smooth com- 
plete intersections of constant multiplicity (see (TS] for similar results). Before 
describing the algorithm, we need a definition which captures this concept. 

Definition 2.15. With the above notation, a dense Zariski-open subscheme U oiS 
such (^^^iU) — >■ U is smooth and free is said to be an /-optimal subscheme of S. 



Corollary 2.16. Let S = and consider the following sequence of instructions. 

(1) Compute D(a,x) = det(Jac^-(a, x)). 

(2) Let J(a, x) = /(a, x) + (-D(a, x)) and compute H — J(a, x) n if [a] . 

(3) Lf H = (0) return "There is no /-smooth subscheme of " and stop. 

(4) Choose h{a) eH\0 and let = A™ \ {a e A™ | h{a) = 0}. 

(5) Choose a term ordering a on T" and compute the reduced a-Grobner basis 
G(a,x) o//(a,x)/f(a)[x] 

(6) Let T = T" \ LTo.(/(a, x)/f (a)[x]), compute the cardinality of T and call 
it fi; then compute the least common multiple of all the denominators of 
the coefficients of the polynomials in G'(a, x), and call it d{a); finally, let 

= A™ \ {a e A™ I d(a) ^ 0} and let U ^UinU2. 

(7) Return Ui, U2, U, T, /i. 

This is an algorithm which returns lAi which is L-smooth, IA2 which is I-free, 11 
which is I-optimal, T which provides a basis as K-vector spaces of all the fibers 
overU2, and fj, which is the multiplicity of all the fibers overU2. 

Proof. It suffices to combine Theorem |2 . 1 2 1 and Proposition |2.6| □ 



Example 2.17. We consider the ideal / — (/i,/2) of K[x,y] where /i = xy — 6, 
f2 = x^+y^ — 13. It is a zero-dimensional complete intersection and we embed it 
into the family /(a, x) = {aixy + 02, a^x^ + O'AV^ + 15)- We compute the reduced 
DegRevLex-Grobner basis of /(a, x)ii'(a)[x] and get 

+ -I- ^, xy^'^, y^ - '^x+ '^y} 

according to the above results, a free locus is given by 010^,04^ 7^ 0. Now we compute 
i?(a, x) = det(JacF(a, x)) and get /^(a, x) = —201030;^ -|- 2a\a/^y'^ . 

We let J(a, x) = i(a, x) -|- (i3(a, x)) and compute J(a, x) n if [a]. We get the 
principal ideal generated by a^a^^a^^ — ^a\a\. In conclusion, an i-optimal sub- 
scheme is W = A^ \ F where F is the closed subscheme defined by the equa- 
tion aiazai[a^azai — \cL\a\) — 0, and /i = 4. 
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Definition 2.18. We say that a point is complex if its coordinates are complex 
numbers, and we say that a point is real if its coordinates are real numbers. 

The following example illustrates the fact that even if we start with a set of 
real points, a zero-dimensional complete intersection which contains them may also 
contain complex non-real points. 

Example 2.19. Let X be the set of the 10 real points {(-1,-1), (2, 8), (-2,-8), 
(3,27), (-3,-27), (4,64), (5,125), (-5,-125), (6,216), (-6,-216)}. A zero-dim- 
ensional complete intersection containing X is {/i,/2} where /i = y — and 
/2 = x2y2 _ i/4095?;4 + 1729/15x2?/ - 74/15x?/2 + l/l5y^ - 8832/5x^ + 5852/15x2/ - 
10754/315y2 + 2160x - 4632/5y -I- 250560/91. Let / denote the vanishing ideal of 
the 10 points and let J denote the ideal generated by {/i, /2}. The colon ideal J : / 
defines the residual intersection. Since J is the intersection of a cubic and a quar- 
tic curve, the residual intersection is a zero-dimensional scheme of multiplicity 2. 
Indeed, a computation (performed with CbCcA) shows that J : / is generated by 
(x -I- 1 /78y - 87/26, y^ - 756y + 658503) . Since 756^ - 4 * 658503 = -2062476 < 0, 
the two extra points on the zero-dimensional complete intersection are complex, 
non real points. 

Theorem 2.20. Let f (x) be a zero- dimensional complete intersection in M[x] and 
let f(a, x) G M[a, x] be a zero- dimensional family containing f(x). Let I be the 
ideal in M[x] generated by f (x), assume that there exists an L-optimal subscheme U 
of , and let aj CzU be the point in the parameter space which corresponds to L. 
If IJ-MJ is the number of distinct real points in the fiber over olj (i.e. zeroes of L), 
then there exist an open semi-algebraic subscheme V oflA such that for every a G V 
the number of real points in the fiber over ol is /ijj./. 

Proof. We consider the ideal I = /(a, x)M(a)[x]. It is zero-dimensional and the 
field ]R(a) is infinite. Since a linear change of coordinates does not change the 
problem, we may assume that I is in x„-normal position (see [T5] , Section 3.7). 



Moreover, we have already observed (see Remark 2.7) that in Proposition |2.6| the 
choice of a is arbitrary. We choose a = Lex and hence the reduced Lex-Grobner 
basis of I has the shape prescribed by the Shape Lemma (see fl5] Theorem 3.7.25). 
Therefore there exists a univariate polynomial hg^ G M(a)[a;„] whose degree is the 
multiplicity of both the generic fiber and the fiber over aj, which is the number 
of complex zeros of /. Due to the shape of the reduced Grobner basis, a point 
is real if and only if its x„-coordinate is real. Therefore it suffices to prove the 
following statement: given a univariate square-free polynomial ft-a S K(a)[xn] such 
that hai has exactly fi^j real roots, there exists an open semi-algebraic subset 
of such that for every point a in it, the polynomial /iq, has exactly /iR^/ real 
roots. This statement follows from [3, Theorem 5.12 where it is shown that for 
every root there exists an open semi-algebraic set in which isolates the root. 
Since complex non- real roots have to occur in conjugate pairs, this implies that real 
roots stay real. □ 

Let us see some examples. 

Example 2.21. We consider the ideal / = {xy — 2y'^ + 2y, x^ — y'^ ~ 2x) in M[x, y], 
and we embed it into the family /(a, x) = (xy — ay^ -\- ay, x^ — _ 2x). We 
compute the reduced Lex-Grobner basis of /(a, x)M(a)[x] and get 

{x^ - 2x - y^, xy - ay^ -f- ay, y^ - ^y^ + ^^^y} 
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Applying the algorithm illustrated in Corollary |2 . 16| we get an /-smooth subscheme 
of for a{a + 2) 7^ 0, and an /-free subscheme for (a — l)(a + 1) 7^ 0. For a 
different from 0,-2, 1,-1 we have an /-optimal subscheme and the multiplicity 
is 4. 

Our ideal / is obtained for a — 2, and hence it lies over the optimal subscheme. 
It has multiplicity 4 and the four zeros are real. 

The computed Lex-Grobner basis does not have the shape prescribed by the 
Shape Lemma, so we perform a linear change of coordinates by setting x — x + 
y, y = X — y. We compute the reduced Lex-Grobner basis and get 

{x-f4^y3_2^^._^y^ y^-y^-\^^y- + \^^y} 
It has the good shape, so we can use the polynomial 

h. = v'~y'-\^,v' + \^,v^ viv - i){y' - \^,) 

We get the following result. 

• For a < —1, a 7^ —2 there are 4 real points. 

• For — 1 < a < there are 2 real points. 

• For a > 0, a 7^ 1 there are 4 real points. 

To complete our analysis, let us see what happens at the bad points 0, —2, 1, —1. 

At the primary decomposition of the ideal /q is {x — 2, j/) n {y^ + 2x,xy,x^), 
hence the fiber consists in the simple point (2,0) and a triple point at (0,0). 

At —2 we see that (a; + 1, y — |)n(a;, j/)n(a; — 2, y^) is the primary decomposition 
of the ideal /_2 7 and hence the fiber consists in the simple point (— |, |), the simple 
point (0,0) and a double point at (2,0). 

At —1 the primary decomposition of the ideal /_i is (x, y)n{x — 2, y), hence the 
fiber consists of the two simple real points (0, 0) and (2, 0). 

At 1 we see that {x,y) H {x — 2,y) Ci {x + \,y — |) is the primary decomposition 
of the ideal Ii, hence the fiber consists of the three simple real points (0, 0), (2, 0), 

V 4' 4^ 

Example 2.22. We consider the ideal / — [xy + 1, -f — 5) in M\x, y], and we 
embed it into the family I[a.,x,y) = {xy + aix + 1, -I- j/^ -I- 02). We compute the 
reduced Lex-Grobner basis of /(a, x)/ir(a)[a;, y] and get G{a,x,y) — {(?i,(?2} where 

gi = X - y^ - aiy'^ - a2y ~ aia2, 

92 = y^ + 2aiy'^ + {a\ + a2)y'^ + 2aia2y + {a\a2 + I) 

which has the shape prescribed by the Shape Lemma (see [15] Theorem 3.7.25). 
There is no condition for the free locus, and D{a.,x,y) = det(JacF(a, x, y)) = 
—2x^ + 2y^ + 2aiy. We let J(a, x, y) = /(a, x,y) -I- {D{a,x,y)) and compute 
J(a, x,y) n K[a\. We get the principal ideal generated by the following poly- 
nomial /i(a) = of a2 -I- Safoj -I- of -I- 3of 03 + 20a^a2 + a| — 8a| -|- 16. An /-optimal 
subscheme is U — Aj| \ F where F is the closed subscheme defined by the equa- 
tion h{a) = 0, and we observe that ^ = 4. 

At this point we know that for /i(a) 7^ each fiber is smooth and has multiplic- 
ity 4, hence it consists of 4 distinct complex points. What about real points? 
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The real curve defined by h{a) = is shown in the above picture. It is the union 
of two branches and the isolated point (0,2). The upper region Rl (with the 
exception of the point (0,2)) corresponds to the ideals in the family whose zeros 
are four complex non-real points. The regions R2 and R3 correspond to the ideals 
whose zeros are two complex non-real points and two real points. The region 
R4 corresponds to the ideals whose zeros are four real points. To describe the 
four regions algebraically, we use the Sturm-Habicht sequence (see [H]) of g2 £ 
M(a)[y]. The leading monomials are y*, 4y^, 4r(a)y^, — 8£(a)?/, 16/i(a) where 
r(a) — al — 2a2, £{a.) = a\a2 + 2a\a\ -I- 2a\ + a\ — 4a2. To get the total number of 
real roots we count the sign changes in the sequence at — oo and +oo; in particular, 
we observe that in the parameter space the ideal I corresponds to the point (0, —5) 
which belongs to the region R4. We get 

R4 {a e I r(a) > 0, l{a.) < 0, h{a) > 0} 



which is semi-algebraic open, not Zariski-open. 



3. Condition Numbers 

In this section we introduce a notion of condition number for zero-dimensional 
smooth complete intersections in ]R[x]; the aim is to give a measure of the sensitivity 
of its real roots with respect to small perturbations of the input data, that is small 
changes of the coefficients of the involved polynomials. 

The section starts with the recall of well-known facts about numerical linear 
algebra. We let m, n be positive integers and let Matnixn(lR) be the set of m x n 
matrices with entries in M; if m = n we simply write Matn(M). 

Definition 3.1. Let M — {rriij) be a matrix in Matmxn(IR)) v = (i;i,...,i;„) a 
vector in M" and || • || a vector norm. 

(a) Let 7' > 1 be a real number; the r-norm on the vector space M" is defined 

by the formula ||t;||.r — (X]"=i biD'' foi' every v G M". 

(b) The infinity norm on M" is defined by the formula ||w|joo = inaxi|tii|. 

(c) The spectral radius q{M) of the matrix M is defined by the formula 
q{M) — maxi I Ail, where the are the complex eigenvalues of A/. 
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(d) The real function defined on Mat,„xn(]R) by M i— > max||^,||=i is a 
matrix norm called the matrix norm induced by || • ||. A matrix norm 
induced by a vector norm is called an induced matrix norm. 

(e) The matrix norm induced by || • ||i is given by the following formula ||Af ||i — 
maxj(^j jmijl). The matrix norm induced by || • ||oo is given by the formula 
||-M||oo = ma'X.j(^j- |TO.y|). Finally, the matrix norm induced by ||-|j2 is given 
by the formula ||Af||2 ~ maxi((Ti) where the Ui are singular values of M. 

If no confusion arises, from now on we will use the symbol || • || to denote both a 
vector norm and a matrix norm. We recall some facts about matrix norms (see for 
instance [1], [T3]). 

Proposition 3.2. Let M he a matrix in Matn(M), let I he the identity matrix of 
type n and let \\ ■ || be an induced matrix norm on Matn(M). // the matrix I + M is 
invertihle then (1 - ||Af||) \\{I + M)-^\\ < 1. 

Proposition 3.3. Let M e Matmxn(II^) and denote by Mi the i-th row of M . Let 
ri > l,r2 > 1 he real numbers such that + = then 

ma.x\\Mi\\r2 < \\M\\r^ < m^/''^ ma.x\\Mi\\r^ 

i i 

In particular, for ri = r2 = 2 

max||Af,;||2 < \\M\\2 < ^/TOmax ||Afi||2 

i i 

This introductory part ends with the recollection of some facts about the poly- 
nomial ring -R'[x]. In particular, given rj — (771, . . . ,77„) S N" we denote by \r]\ the 
number r]i + . . . + rjn, by 77! the number 771! . . . and by x** the power product 

■^1 ■ ■ ■ -^n ■ 

Definition 3.4. Let p be a point of X"; the if-linear map on i^[x] defined by 
/ ^ f{p) is called the evaluation map associated to p and denoted by evp{f). 

Definition 3.5. Let c? be a nonnegative integer, let r > 1 be a real number, let p 
be a point of M" and let ^(x) be a polynomial in M[x]. 

(a) The formal Taylor expansion of g(x) at p is given by the following expres- 
sion: 5W = E|,,i>o 

(b) The polynomial E|,,|>d - pY' is denoted by g^'^{x,p). 

(c) The r-norm of ^(x) at p is defined as the r-norm of the vector |f (p). If 

(p)lli- = 1 then g{x) is called unitary at p. 
We use the following formulation of Taylor's theorem. 

Proposition 3.6. Let p be a point o/M" and let ^(x) be a polynomial in IR[x]. For 

every point q G M" we have 

giq) ^ gip) + Jacg(p)(g -p) + ^{q~ pYHg{0{q - p) 

where is a point of the line connecting p to q and Hg{^) is the Hessian matrix 
of g at 

Given f (x) = {/i(x), . . . , /„(x)}, a zero-dimensional smooth complete intersec- 
tion in K[x], we introduce a notion of admissible perturbation of f(x). Roughly 
speaking, the polynomial set £(x) — {ei(x), . . . , e„(x)} C M[x] is considered to 
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be an admissible perturbation of f(x) if the real solutions of (f + £)(x) = are 
nonsingular and derive from perturbations of the real solutions of f(x) — 0. Using 
the results of Section [2] we formalize this concept as follows. 

Definition 3.7. Let f (x) = {/i(x), . . . , /„(x)} be a zero-dimensional smooth com- 
plete intersection in M[x], let /iRj be the number of real solutions of f (x) — 0, and 
let £(x) = {ei(x), . . . ,e„(x)} be a set of polynomials in M[x]. Suppose that the 
assumptions of Theorem 2.20 are satisfied, let V C A]^ be an open semi-algebraic 
subset of U such that a/ e V, and for every a e V the number of real roots of 
f(Q;,x) = is equal to /iK,/. If there exists a e V such that (f + £)(x) ~ f(a,x), 
then e(x) is called an admissible perturbation of f(x) . 

Henceforth we let £(x) — {ei(x), . . . , £„(x)} be an admissible perturbation 
of f(x), and let Z^(X) = {pi, . . . ZR(f + e) = {ri, . . . , r^g ,} be the sets 

of real solutions of f(x) = and (f + £)(x) = respectively. We consider each ri 
as a perturbation of the root p^, hence we write — pi + Apj for i = 1, . . . , fi^j- 

Now we concentrate on a single element p of ^^(f). 

Corollary 3.8. Let p be one of the real solutions off = 0, and p + Ap the corre- 
sponding real solution of t + e — 0. The we have 

(1) 0= (f + £)(p + Ap) =£(p) + Jacf+,(p)Ap+(t-i(a),...,t-„(^„))* 

where ^i, . . . , ^„ are points on the line which connects the points p and p + Ap, and 
"^ji^j) = l^P^Hf^+^. {£,j)Ap for each j = 1, . . . , n. 



Proof. It suffices to put q = p + Ap, apply the formula of Proposition 3.6 to the 
polynomial system (f -I- £)(x), and use the fact that f (p) = 0. □ 

Example 3.9. We consider the zero-dimensional smooth complete intersection 
f = {/i,/2} where /i — xy — 6, f2 = x'^ + y^ — 13 and observe that ZR(f) = 
{(-3, -2), (3, 2), (-2, -3), (2, 3)}. The set f (x) is embedded into the following fam- 
ily F(a,x) = {xy + ai,x'^ + a2y^ + a^}. 
The semi-algebraic open set 

V = {aeW^\al- \a\a2 > 0, a2 > 0, as < 0} 

is a subset of the /-optimal scheme U = {oi ct2(oi\ — 4af 02) 7^ 0}. Moreover, 

it contains the point ct/ — (—6, 1,-13), and the fiber over each a G V consists 
of 4 real points. The set £(x) — {S\,82y^ 4- (Ja}, with bi E R, is an admissible 
perturbation of f (x) if and only if the conditions {S3 — 13)^ — 4((5i — 6)^(62 + 1) > 0, 
$2 > —1, and 63 < 13 are satisfied. Since the values 5i — 2, S2 — ^, and (53 = 

satisfy the previous conditions, the polynomial set £(x) = {2, |?/^} is an admissible 
perturbation of f (x). The real roots of (f + £)(x) = are 

Z«(f + £) = { (-3, -I) , (3, I) , (-2, -2), (2, 2)} 

For each <E -ZR(f + £) the matrix Jacf+e(ri) is invertible, as predicted by the 
theory. On the contrary, by evaluating JaCf+e(x) at the third and the fourth point 
of 2'M.{f) we obtain a singular matrix. This is an obstruction to the development of 
the theory which suggests further restrictions (see the following discussion). 



Our idea is to evaluate Ap using equation ([T]) of Corollary 3.8 However, while 



the assumption that £(x) is an admissible perturbation of f(x) combined with the 
Jacobian criterion guarantee the non singularity of the matrix JaCf+£(p -I- Ap), 
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they do not imply the non smgularity of the matrix Jacf+£(p), as we have just 
seen in Example ] 3. 9 1 The next step is to find a criterion which guarantees the non 
singularity of Jacf+£(p). 

Lemma 3.10. // || Jacf (p)"-'^ Jac£(p)|| < 1 then Jacf_|_£(p) is invertible. 

Proof. By assumption p is a nonsingular root of f(x) = 0, hence JaCf(p) is in- 
vertible and so JaCf+E:(p) can be rewritten as JaCf+e(p) = JaCf(p) + JaCe(p) = 
Jacf{p) (/ + Jacf (p)^^ JaCe(p)). Consequently, it suffices to show that the matrix 
/ + JaCf JaCe(p) is invertible. And we achieve it by proving that the spec- 
tral radius e(JaCf JaCe(p)) is smaller than 1. We have e(JaCf JaCe(p)) < 
II JaCf Jac£(p)|| < f , and the proof is now complete. □ 

Note that the requirement || JaCf (p)~^ JaCe(p)|| < f gives a restriction on the 
admissible choices of £(x), as we see in the following example. 



Example 3.11. (Example 3.9 continued) 

Let £(x) ~ {Si,S2y^ + S^}, with 6i e M, be an admissible perturbation of the 



zero-dimensional complete intersection f (x) of Example 3.9 We consider the real 
crvlntirvn n , — ^"^,3) of f = and computc II Jacf (p4)^^ Jac£(p4)||| ~ ly'^l- From 
the condition \62\ < ^y/T3 is sufficient to have JaCf+£(p4) invertible. 



Lemma 



3.f0 



From now on we assume that the hypothesis of Lemma 3.10 is satisfied. In order 
to deduce an upper bound for || Ap|| we consider an approximation of it. 

Definition 3.12. If || JaCf (p)^^ JacE(p)|| is different from 1, we denote the number 
1/(1 — II Jacf JaCe(p)ll) by A(f, e,p). Moreover, if equation (fl]) is truncated 
at the first order, we get the approximate solution ~ 3acf^s{p)^^{p) which we 
call Ap^ . 

Proposition 3.13. Assume that \\ JaCf (p)^"'^ JaCe(p)|| < 1 and let \\ ■ \\ be an induced 
matrix norm. Then we have 

(2) ||Api|| <A(f,£,p) ||Jacf(p)-i|| ||£(p)|| 



Proof. Lemma 3.10 guarantees that the matrix Jacf^^^p) is invertible, so 
Ap^ = - Jacf+e(p)"^£(p) = -(Jacf (p) + Ja.Cs{p)r^e{p) 
= - (/ + Jacf(p)"^ Jac£(p)) ^ Jacf(p)~^£(p) 



We apply the inequality of Proposition 3.2 to Jacf (p) Ja.Cs{p), and get 



||Api|| < ||(/-fJacf(p)'iJac,(p)-i)|| ||Jacf(p)-i|| ||£(p)|| 
< A(f,£,p) ||Jacf(p)-i|| ||£(p)|| 

which concludes the proof. □ 

We introduce the local condition number of the polynomial system f (x) = 0. 

Definition 3.14. Let f(x) be a zero-dimensional smooth complete intersection 
in M[x], let p be a nonsingular real solution of f (x) = 0, and let || • || be a norm. 

(a) The number K{i,p) = \\ Jacf (p)~^|| || Jacf(p)|| is called the local condition 
number of f (x) at p. 

(b) If the norm is an r-norm, the local condition number is denoted by Kr{i,p). 
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The following theorem illustrates the importance of the local condition number. 
It depends on / and p, not on e and is a key ingredient to provide an upper bound 

for the relative error 



Theorem 3.15. (Local Condition Number) 

Let II • II be an induced matrix norm; under the above assumptions and the condition 
II Jacf (p)"-"^ JaCe(p)|| < 1 we have 

(3) ^<A(f,s,,).(f,p)^l|J^^^(^)ll . 11^(0) -^^(O'^)ll 



P\\ - ' ' ' ' ' Vl|JaCf(p)|| ||f(0)-f>2(0,p)|| 



Proof. By Definition |3.5| the evaluation of e at can be expressed in this way 
£(0) = £{p) - 3aCe{p)p + £-^(0,p), and so £{p) — £(0) + 3aCs{p)p — £-'^{0,p). 
Dividing ^ of Proposition 3.13 by ||p|| we obtain 

ll^^^ll < A(f,£,,)||JacKp)-"^(^)" 



IN - ^ ' " IHI 

JaCe(p)|||b|| + ||£(0)-£^2(o^p)|| 



< A(f,£,p) ||Jacf(p)-i| 



A(f,£,p)|| Jacf(p)-i|| ||Jac,(p)| 



llPlI 

|g(0)_e>2(o,p)|| 

Ibl! 



Using again Definition 3.5 we express f(0) = i{p) — Jacf(p)p + f-^(0,p); since 
f(p) we have ||f(0) — f-^(0,p)|| = || Jacf(p)p|| < || Jacf (p)|| ||p|| from which 

1 ^ ||JaCf(p)|| 



IIpII - ||f(0)-f>2(0,p)|| 
We combine the inequalities to obtain 

< A(f,£,p)|| Jacf(p)-^|| (ll Jac,(p)|| + || Jacfb)|| ~ ^^l^^^^^^ 

Ml. /X Inn. / ^n /II JaCe(p)|| 1 1 £(0) - E^^ (0, ja) || 

< A(f,£,p)|| Jacf(p)-i|||| Jacf(p)r " ■ " ^ 



Jacf(p)|| ||f(0)-f>2(0,p)|| 
and the proof is concluded. □ 

The following remark contains observations about the local condition number. 

Remark 3.16. We call attention to the following observations. 

(a) The notion of local condition number given in Definition 3.14 is a generaliza- 
tion of the classical notion of condition number of linear systems (see jj). In 
fact, if f (x) is linear, that is f (x) = Ax—b with A e Matn(M) invertible, and 
Z-R{f) = {p} = {A^^b}, then K{f,p) is the classical condition number of the 
matrix A. In fact JaCf(x) = A, and so «;(f,_p) = || JaCf (p)~^|| || JaCf(p)|| = 
||yl^^||||A||. Further, if we consider the perturbation £(x) = AAx — Ab, 
relation (|3| becomes 

< ' ILl-^ll IIAII ('"^^" + "^''"^ (4) 

which is the relation that quantifies the sensitivity of the Ax = b problem 
(see [1], Theorem 4.1). 
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(b) Using any induced matrix norm, the condition number K(f,p) turns out 
to be greater than or equal to 1. In particular, using the 2-norm we have 
K2{{,p) — j"'""'(jacf (pjj ' ^^^^ local condition number attains its 
minimum, that is K2(fj?') — Ij when JaCf(p) is orthonormal. 

(c) The condition number K{i,p) is invariant under a scalar multiplication 
of the polynomial system f(x) by a unique nonzero real number 7. On 
the contrary, K(t,p) is not invariant under a generic scalar multiplica- 
tion of each polynomial /j (x) of f(x). The reason is that if we multiply 
each /j(x) by a nonzero real number jj we obtain the new polynomial set 
g(x) = {71/1 (x), . . . ,7„/„(x)} whose condition number at p is 

K{g,p) - II Jacf(p)-ir-i|| ||r Jacf(p)|| ^ k{{,p) 

where F — diag(7i, . . . , 7„) g Matn(M) is the diagonal matrix with en- 
tries 71, . . . ,7„. 

(d) It is interesting to observe that if p is the origin then Formula (|3| of the 
theorem is not applicable. However, one can translate p away from the ori- 
gin, and the nice thing is that the local condition number does not change. 



4. Optimization of the local condition number 

In this section we introduce a strategy to improve the numerical stability of 
zero-dimensional smooth complete intersections. Let f (x) = {/i(x), . . . , /„(x)} be 
a zero-dimensional smooth complete intersection in K[x], and let / be the ideal 
of M[x] generated by f (x); our aim is to find an alternative representation of / with 
minimal local condition number. 



Motivated by Remark 3.16 item (b) and (c), we consider the strategy of resiz- 



ing each polynomial of f(x), and study its effects on the condition number. The 
following proposition shows that rescaling each /j (x) so that j^{p) has unitary 
norm is a nearly optimal, in some cases optimal, strategy. The result is obtained 
by adapting the method of Van der Sluis (see [TSj, Section 7.3) to the polynomial 
case. 

Proposition 4.1. Letp he a nonsingular real solution o/f (x) = 0, let ri > 1, r2 > 1 

he real numbers such that ^ -|- ^ = 1, including the pairs (l,c») and (cxd, 1), let 
7 — (71,..., 7„) he an n-tuple of nonzero real numhers, and let g^(x), m(x) he 
the polynomial systems defined by g'y(x) — {7i/i(x), . . . , 7„/„(x)} and it(x) = 

{ll|^(p)L-//lW,---Jlllf(P)llr//nW}- 

(a) We have the inequality Kr-^{u,p) < n^^^^ Kr-^{g'y,p) . 

(b) In particular, if {ri,r2) = (00, 1) we have the equality 

Kooiu,p) = min.^Koo(g7,p) 
where «(x) = {|| |^(p)||- ^^(x), . . . , \\^{p)\\^^ f,,{^)}. 

Proof Let F = diag(7i, . . . , 7„) and D = diag(|| |f^(p)||-i, . . . , |||^(p)||-i); then 
JaCg (x) — F JaCf(x) and JaCit(x) = Jacf(x). The condition numbers of g-,.(x) 
and it(x) at p are given by 

p) = ||(FJaCf(p)) i||,rJ|F Jacf(p)||ri 

Kr,iu,p) = \\{DJa.Cf{p)y^\\r^\\DJaCf{p)\\r, 
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From Proposition |3.3| we have 



\\D 3acf{p)\\r, 

\{D3SiCf{p)r^\\r, 



< 



< 



JaCf 



\{D 3a.Cf{p)),\ 



f \P)D- 



Jac, 



9x 



ip) 



< ||Jac^\p)r" 



r Jacf(p)||, 



therefore Kr^{u,p) < (g-,.,p) and (a) is proved. To prove (5) it suffices to 

use (a) and observe that v}/^ = 1 □ 

Remark 4.2. The above proposition imphes that the strategy of rescahng each 



polynomial fj (x) to make it unitary at p (see Definition 3.5 1 is beneficial for lowering 
the local condition number of f (x) at p. This number is minimum when r = oo, 
it is within factor yjn of the minimum when r = 2. However, for r = 2 we can 
do better, at least when all the polynomials /i(x), . . . , /„(x) have equal degree. 



The idea is to use Remark 3.16 item (b) which says that when using the matrix 



2-norm, the local condition number attains its minimum when the Jacobian matrix 
is orthonormal. 

Proposition 4.3. Let f = (/i,...,/„) be a smooth zero- dimensional complete 
intersection mM[x] such that deg{fi) — ■ ■ ■ — deg(/„) and letp G Z^{i). Moreover, 
let C = (cij) G Mat„(M) be an invertible matrix, and let g be defined by g'"^ — C-f*"^. 
Then the following conditions are equivalent 

(a) K2(giP) = 1; the minimum possible. 

(b) C*C = (Jacf(p) Jacf 

Proof. We know that K2(g,p) = 1 if and only if the matrix JaCg(p) is orthonormal. 
This condition can be expressed by the equality JaCg(p) JaCg(p)* = /„, that is 
C Jacf (p) Jacf (p)*C* — /„ and the conclusion follows. □ 

We observe that condition (b) of the proposition requires that the entries of C 
satisfy an underdetermined system of {n'^ + n)/2 independent quadratic equations 
in unknowns. 



5. Experiments 



In numerical linear algebra it is well-known (see for instance [4], Ch. 4, Section 1) 
that the upper bound given by the classical formula (4) of Remark 3.16 (a) is not 
necessarily sharp. Since our upper bound ([s]) generalizes the classical one, as shown 
in Remark |3.16[ we provide some experimental evidence that lowering the condition 
number not only sharpens the upper bound, but indeed stabilizes the solution point. 
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Example 5.1. We consider the ideal / = (/i, /2) in 

/i = |a;2y + a;?;2 + iy3 + ia;2- 

h ~ X + j^xy + ^x — 52^2/ I i3y 52 

It is a zero-dimensional smooth complete intersection with 7 real roots and we 
consider the point p = (0,1) G -ZR(f). The polynomial system f ~ {/i,/2} is 
unitary at p and its condition number is K2(f,p) 



_ 4_ 
13 



8. Using Proposition 4.3 



we 



construct a new polynomial system g with minimal local condition number at p. 
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The new pair of generators g is defined (see Proposition 4.3 ) by the following the 
formula g*'' = C • f*'', where C = (cy) G Mat2(M) is an invertible matrix whose 
entries satisfy the following system 



^2 , „2 -25 

"--11 + "^21 — 16 

C11C12 + C21C22 = -J 

'-12 + "^22 — 16 



A solution is given by Cn = 1, C12 = 0, C21 — f|, C22 — ^ff, and we observe 
that the associated unitary polynomial system g = {/i, f|/i — f§/2} provides an 
alternative representation of / with minimal local condition number K2{g,p) = 1 
at the point p. 

Now we embed the system f (x, y) into the family F{a, x, y) — {Fi, ^"2} where 
Fi{a,x,y) = \x''y + xy'' + \y^ + lx'' -lxy+{f^~a)y'' 

+ {l+a)x+{-l+a)y+l^-2a 
F2{a,x,y) = + ^xy^ + ~ §xy + {f^ + a) y^ 



+ {-^ + a) x - + a) y ~ + a 



J_ I „2 
13 



We denote by /^(a, x, y) the ideal generated by F{a, x, y) in M[a, a;, y], compute the 
reduced Lex-Grobner basis of If{cl, x, y)^{a)[x, y\, and get 

where li{a,y),l2{a,y) € K[a,y] have degree 8 in y and dp{a) € M[a] has degree 12. 
This basis has the shape prescribed by the Shape Lemma and a flat locus is given 
by {a S K I dpia) ^ 0}. We let Dp{a,x,y) = det(Jac_F(a, a;, y)), JF{a,x,y) = 
/^(a, X, y) + {Dp {a, x,y)), compute Jpia, x,y)r]M.[a], and we get the principal ideal 
generated by a univariate polynomial hF{a) of degree 28. An /-optimal subscheme 
is Up = {a G R I dp{a)hp{a) 7^ 0}. An open semi-algebraic subset Vp oiUp which 
contains the point ai — {) and such that the fiber over each a ^ Vp consists of 7 
real points, is given by the open interval (ai, 0.2), where ai < and 02 > are the 
real roots of dp{a)hp{a) — closest to the origin. Their approximate values are 
ai = -0.00006 and 02 = 0.01136. 

To produce similar perturbations, we embed the system g{x, y) into the family 
G{a,x,y) = {G'i,G2} where 

Gi{a,x,y) = lx^y + xy^ + \y^ + lx^-lxy+{§-a)y^ 

+ {i^+a)x+{~l + a)y+i^-2a 
^2(0, X, y) = ~Te,'^ "t" 64"^ y ~ T6''^y ~^ eiV sW'^ ~ 128'^^ 

+ i-m + -)y' + i'ii^ + -)- + i-fo--)y + m + -' 

We denote by Ida, x, y) the ideal generated by G{a, x, y) in IR[a, x, y], compute the 
reduced Lex-Grobner basis of /G(a, x, y)M(a)[a;, y], and get 

{^ + ^' y' + h{a,y)} 

where l3{a,y),l4{a,y) G M[a,?/] have degree 8 in y and dc{a) G R[a] has degree 12, 
therefore the basis has the shape prescribed by the Shape Lemma. A flat lo- 
cus is given by {a G M | rfcC") 7^ 0}. We let DG{a,x,y) = det(JacG(a, a;, 2/)), 
JG[a,x,y) — lG{a,x,y) + {DQ{a,x,y)) and compute JGia,x,y) n ]R[a]. We get 
the principal ideal generated by a univariate polynomial hcio) of degree 28. An 
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/-optimal subscheme is Ug = {a G M | dG{a)hG{a) 7^ 0}. An open semi-algebraic 
subset Vg of Ug containing the point ai — i) and such that the fiber over each 
a G Vg consists of 7 real points is given by the open interval (as, 04), where 03 < 
and a4 > are the real roots of dG{a)hG{a) = closest to the origin. Their 
approximate values are 0:3 = —0.00009 and Q!4 ~ 0.00914. 

(ai,a4). According to Definition 3.7 the polynomial set £{x,y) — 



Let a G 
,,2 I 



{—ay^ + ax + ay — 2a, ay^ -\- ax — ay + a'^} is an admissible perturbation 



of f{x,y) and g{x,y). Further, since || Ja cf (p) ^JaCe(p)||2 = -v/BSlaj < 1 and 



II JaCg(p) ^ JaCe(p)||2 = V^\ct\ < 1 Theorem 3.15 can be applied 



We let q S Zr{{ + e) and r £ Zjg^g + e) be the two perturbations of the point p. 
In order to compare the numerical behaviour of f and g at the real root p we 
compare the relative errors and ^^Tr^^^" for different values of a. The first 

^ \\P\\2 \\P\\2 

column of the following table contains the values of the local condition numbers 
of f and g at p. The second column contains the mean values of the upper bounds 



UB(f,p) and UB(g,p) given by Theorem 3.15 computed for 100 random values 
of a € (ai,a4). The third column contains the mean values of and 
for the same values of a. 



K2(f,p) 


UB(f,p) 




8 


0.1729 


0.000097 


l^2{g,p) 


UB(g,p) 


II' — PII2 
IIpIU 


1 


0.0275 


0.000023 



The fact that the mean values of ^^'?| i^^^'' are smaller than the mean values of u^^^ 

\\ph Wvh 

suggests that p is more stable when it is considered as a root of g instead of as a 
root of f . 



17-^ 
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xy 

70, 



ill, f2, fs 
24 



in M[a;, y, z] where 



89^2^ + 
^ + 2xz + 



tX 



yz 



39t 
89-^ 



_6_ 
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f 



Example 5.2. We consider the ideal / 
/i 
h 

fs = y 

It is a zero-dimensional smooth complete intersection with 6 real roots and we 
consider the point p — (1, 0, 0) £ 2k (f). The polynomial system f = {/i, /2, /s} is 
unitary at p and its condition number is K2(f,p) — 123. Using Proposition |4.3| we 
construct a new polynomial system g with minimal local condition number at p. 
The new set g is defined by g*'' = C - f'"^, where C = {cij) £ Mat3(IR) is an invertible 
matrix whose entries satisfy the following system 



--11 



c: 



21 



c 



31 



C11C12 -I- C21C22 
C11C13 -I- C21C23 

"^12 ^ ''22 ^ ''32 
'3l2'^13 + C22C23 
'^13 + '^23 + C33 



A solution is given by cn 



'333 



C31C32 
C31C33 

'332'333 
1, C12 



57229225 

ileo 

15129 



57229225 
15129 




Cl3 



C23 



'^32 



0, C21 



C22 



7565 
123 ■ 



7564 
123 ' 



Therefore the associated unitary polynomial system is the following 
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g — {/i, ^Y^fi — '^f§'f2, fa}- It provides an alternative representation of / with 
minimal local condition number K2{g,p) = 1 at the point p. 

We embed the system f (a;, y, z) into the family F{a, x, y, z) = {Fi, F2, F^} where 

Fi{a,x,y,z) = ^x'' + (l ~ a^)xy + + a)x + {-^ - a)y + {-^ + a^) 
F2ia,x,y,z) = + + a)xy + yz + {§ + a)x + + a)y 

F3{a,x,y,z) = y'^ + 2xz + {1 - 2a)z'^ + {-1 + a)z 

We denote by If{o-, 2;, y, z) the ideal generated by -F(a, x, y, z) in M[a, x, y, z], com- 
pute the reduced Lex-Grobner basis of If{ci, x, y, z)]R(a)[x, y, z], and get 

l-^ ^ dF(a) ' y dpia) ' ^ 6^(0) J 

where /i(a, z), ^2(0, z), Z3(a, z) G M[a, z] have degrees deg2(/i) = deg^(^2) — 7 and 
deg^(Z3) = 8 while dpia) e M[a] has degree 54, and e_F(a) G M[a] has degree 11. 
The basis has the shape prescribed by the Shape Lemma. A flat locus is given 
by {a e M I dF(oc)e.F{(^) 7^ 0}. We let DF{a,x,y,z) = det(JacF(a, a:, y, 2)), 
7^(0, X, y, z) — iFio-, x,y, z) + {DF^a, x,y, z)) and compute JF(a, x, y, z)nM[a]. We 
get the principal ideal generated by a univariate polynomial hp (a) of degree 59. 
An /-optimal subscheme is Uf = {a € K | dF{a)eF{a)hF{ct) ^ 0}. An open semi- 
algebraic subset Vf of Up containing the point a/ = and such that the fiber over 
each a G Vf consists of 6 real points is given by the open interval (ai, a2), where 
«! < and a2 > are the real roots of dF{a)eF{a)hF{a) = closest to the origin. 
Their approximate values are ai = —0.17082 and a2 = 0.20711. 

To produce similar perturbations, we embed the system g(a;, y, z) into the family 
G{a,x,y,z) = {Gi, 02,03} where 

Oiia,x,y) ^ ^x'' + {1- a^)xy+i-§+a)x + {-§.- a)y + {^§.+a^) 

02{a,x,y) = _||^a;2 + (5M+a)x2;-I|;^yz + (f|f +a)x 

44254 ,„\y_ 15128 
3485 ^ 3485 

03{a,x,y) = y^ + 2xz + {l-2a)z^ + {-l + a)z 

We denote by Icio-, x, y, z) the ideal generated by G(a, x, y, z) in M[a, x, y, z], com- 
pute the reduced Lex-Grobner basis of Ig{o-, y, z)]R(a)[x, y, z], and get 

I '4(g,z) „, I '5(g,z) „9 I Wo^i 
dG(a) ' y + dcCa) ' ^ + eG(a) J 

where li{a, z),l^(a, z),lQ{a, z) € M[a, 2] have degrees deg2(/4) = deg^(^5) = 7 and 
deg^{l%) = 8 while daia) & K[a] has degree 54, and 6^(0) € K[a] has degree 11. 
The basis has the shape prescribed by the Shape Lemma. A flat locus is given 
by {a £ K I dGi{a)dG2{ct) 7^ 0}. We let DG{a,x,y, z) — det(JacG(a, a;, y, z)), 
Jg{o-, x,y, z) — Ig{0', x,y,z) + {DG{a, x,y,z)) and compute JG{a,x,y, z) nM[a]. We 
get the principal ideal generated by a univariate polynomial hG{<i) of degree 59. 
An /-optimal subscheme is Ug = {a S M | dG{a)eG[oi)hG{a) ^ 0}. An open semi- 
algebraic subset Vg of Ug containing the point a/ = and such that the fiber over 
each a £ Vg consists of 6 real points is given by the open interval (013, 04), where 
a3 < and a4 > are the real roots of dG{a)eG{o)hG{a) = closest to the origin. 
Their approximate values are — —0.02942 and = 0.03312. 



Let a € {a^.oiiC). According to Definition 3.7 the polynomial set £{x,y) = 
{—a^xy + ax — ay + a'^, axy + ax + ay, —2az'^+az} is an admissible perturbation 
of f(a;, y,z) and g(x,y, z). 
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We let q g 2^R(f + e) and r € ZR(g + e) be the two perturbations of the point p. 

In order to compare the numerical behaviour of f and g at the real root p we 

compare the relative errors ^^li"!!!!^^ and ^^jj^y'^^^ for different values of a. The first 

column of the following table contains the values of the local condition numbers 

of f and K at p. The second column contains the mean values of and ^^Tr,?^^^ 

^ llPlb ||p||2 

for 100 random values of a G (01,04). 



«2(f,p) 


h-ph 


Mb 


123 


0.0436 


K2(g,P) 


I''"-P||2 


IIpIU 


1 


0.0221 



As in the example before, the fact that the mean values of i, ,1 ^ are smaller than 

the mean values of ^^'|'|~y'|^^ suggests that p is more stable when it is considered as a 
root of g instead of as a root of f . 
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